Dimension dependence of the conductance distribution in the non-metallic regimes 
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The conductance distribution P(g) in disordered three- 
dimensional systems in the critical and insulating regimes is 
calculated numerically and compared with P(g) of weakly dis- 
ordered quasi-one dimensional systems. Although both sys- 
tem exhibit similar transport properties, their conductance 
distributions differ quantitatively. We explain this difference 
by analysis of the spectrum of transfer matrix. Our numerical 
data confirm that the statistics of the conductance is consis- 
tent with the one-parameter scaling also in the strongly lo- 
calized regime. 

PACS numbers: 71.30. +h, 7I.23.-k, 72.15.Rn 



I. INTRODUCTION 

Our understanding of the metal - insulator transition 
(MIT) in disordered solids is based on the scaling theory 
of localizationEl. We believe that MIT is universal. It 
means that the system-size dependence of any quantity, 
say conductance <?, is determined only by one parameter. 
The validity of the one parameter scaling was tested and 
confirmed in the neighborhood of the critical point by 
various numerical simulations on quasi-one dimensional 
(Q1D) systemsa'0. However, a detailed analytical the- 
ory of MIT is still missing. Their formulation is compli- 
cated also by absence of the self-averaging of the conduc- 
tance in disordered systems at zero temperatureu. Com- 
plete description of transport properties requires there- 
fore a detailed analysis of the conductance distribution 
P{g)- Such theory exists only.in the weak disorder limit. 
The weak disorder expansion!], Dorokhov-Mello-Pereyra- 
Kumar equation (DMPK)Q, and the random matrix the- 
ory (RMT)lIu provide us with entire quantitative descrip- 
tion of the transport statistics. However, critical regime 
never appears in the weak disorder limit. 

For the critical regime, analytical results for the 
conductance statistics were derived only in dimensions 
slightly, above the lower critical dimension d = 2 + e with 
e <c lBE2l. These cannot be. extrapolated to three di- 
mensional systems (e = 1)EH. Our information about 
the critical regime in 3D systems is therefore based only 
on numerical simulations. The shape of the conductance 
distribution at the critical point was analyzedOtJ E3 and 
one-parameter scaling of the mean conductance was nu- 
merically proven inO. 

Recently, an attempt has been made to obtain rele- 
vant information about the critical regiese-from studies 
of the weakly disordered Q1D systemsc3'El. Although 
Q1D systems do not exhibit critical phenomena (corre- 



lation length never diverges ) their transport properties 
resemble that observed in 3D systems: (1) in the metal- 
lic regime when the system length L z is shorter than 
the localization length £, both RMT and the DMPK 
equation predict that P(g) in Q1D systems is Gaussian 
with a constant, disorder independent width. This was 
nn mpri pally observed also for squares (2D)Ej and cubes 
(3D)&E3liil. Onkijthe variance of conductance depends 
on the dimensionB. (2) In the limit of very long system 
length, L z ^> £, Q1D samples become insulating, in spite 
of weak disorder. RMT predicts that the distribution of 
log g is Gaussian and var log g = —2 (log g)u. Log-normal 
distribution wasrobserved also in 2D and 3D strongly dis- 
ordered systemsc3'E3. The dimension-dependence of the 
variance var log g is, however, not understood. (3) In the 
intermediate regime, where (g) ~ 1, analytical, relation 
for the conductance distribution can be obtainedoO. Its 
shape is similar talhat obtained numerically at the crit- 
ical point of MITti 

In this paper we discuss the differences between the 
conductance distribution for Q1D weakly disordered 
systems and for 3D cubes. As the metallicr-jlimit is 
already goad .-understood both analyticallyB~BE3 and 
numericallyE3E3 we concentrate to the critical regime 
((g) ~ 1) and the strongly localized regime ((g) << 1. 
Although the distributions are qualitatively similar in 
each transport regime, quantitative differences will be 
found. We show that statistics of the conductance de- 
pends on the dimension in the critical and even in the 
localized regime. This must be taken into account in 
the entire description of the transport. Our analysis of 
the statistical properties of parameters z, which are de- 
fined by Eq. (0) below, confirms the validity of the one- 
parameter scaling theory in the insulating regime. 

In Section 2 we describe the model we use. Section 3 
presents numerical data for the conductance distribution 
in the intermediate (critical) regime. The shape of the 
distribution in the three dimensional localized regime is 
studied in detail in Section 4. 



II. THE MODEL 

We consider Anderson model defined by the Hamilto- 
nian 

H = Wj2 £ n\n){n\ + Y, \n)W\ + \n')(n\ (1) 



where n counts the sites in d-dimensional lattice and nn' 
are nearest- neighbor sites. The size of the lattice is either 



L d (d — 2,3) or L d ~ x x L z . Random energies s n are 
distributed with the box distribution \e n \ < 1/2 and W 
is the strength of the disorder. 

Our analysis of the conductance distribution is based 
on the Landauer formula for the conductances 



N 



g = Tr th = Y^ cosh _2 (z l /2) 



(2) 



In (0), t is the transmission matrix, and the parameters 
Zi, i. = 1, 2, . . . , N, the parameterize the eigenvalues of 
tH 13. In the limit L z >> L, parameters z converge to 
Lyapunov exponents of the transfer matrix. TV is the 
number of open channelsc^l. 

According to (g) , statistics of the conductance is com- 
pletely determined by statistical properties of the param- 
eters z. It is well known that in weakly disordered sys- 
tems the spectrum of z is linear 

zi cui (3) 

independently on dimension and shape of the sampletMl 
The linear behaviopiH) can be derived alscvirom the ran- 
dom matrix model 2 .] and was used in Ref.c3 to prove of 
the universality of conductance fluctuation^. 

It is important to mention that weakness of the disor- 
der is not only necessary but also sufficient condition for 
the linear form of the spectra (£5|). Relation (Q) holds 
therefore for any weakly disordered-jXlID systems in- 
dependently on the system lengtha'Lr&Ej. This is not 
true when disorder increases. Thus, we have numerically 
proven that 



„d-i 



W = W r 



(4) 



at the critical point of MIT in_dimensiop_(i = 3 and 4. Re- 
lation (|) holds for both Q1DE2 and 3DB systems. In 3D, 
it enables us to understand the dimension-dependence of 
the conductance distribution at the critical pointESEj. 

The form of the spectra of z for strongly disordered 
systems is known for the Q1D samples^. An attempt 
was madcJx) study z for 2D and 3D strongly disordered 
systemsE§c3. We present here much more accurate data 
for parameters z, which enable us to estimate the contri- 
bution of higher channels to the conductance. 



III. CONDUCTANCE DISTRIBUTION IN THE 
INTERMEDIATE (CRITICAL) REGIME 

Three-dimensional Anderson model with a box dis- 
tributed disorder exhibits MIT at the critical point 
W c = 16.5. We calculated the conductance distribution 
for ensembles of three-dimensional cubes with disorder 
10.6 < W < 21.5. The system size was 8x8x8 which 
is large enough to exhibit quantitative properties of the 
conductance distribution. Calculated shape of P(g) is 
compared with the conductance distribution for weakly 



disordered Q1D systems of the size 8 2 x L z . The disorder 
W = 4 assures the existence of the metallic regime for 
short L z , L z < 100. The localized regime appears when 
L z exceeds 250. Tuning the length L z we find for each 
3D ensemble the corresponding Q1D ensemble with the 
same mean conductance. 

A comparison of the conductance distribution of 3D 
and Q1D systems is presented in Figures [y, g and [|. 
In the metallic regime, when (g) > 1, P(g) is Gaussian 
in both 3D and Q1D geometry, as supposed (Fig. pi). 
The width of the distribution is universal, independent 
on the system size. It depends only on the system shape, 
and is slightly smaller for Q1D th&H-,for 3D systems, in 
agreement with theoretical results&EZl. 

Figures 1 (b-e) present the shape of P{g) in the critical 
regime, when (g) m 1. We know that P(g) for 3D systems 
becomes non-anal ytica l at g = 1 and decreases exponen- 
tially when g > lEBIij. The same is-J-rue also for Q1D 
systems in the intermediate regimcEErcd. P(g) decreases 
for g > 1 much faster than predicted theoreticallyc£l. The 
exponential decrease in Q1D is even faster than that in 
3D systems. This is because for a given (zi) the dif- 
ference (22 — z\) is much smaller in 3D than in Q1D, 
as follows from the comparison of (0) with (0). Higher 
channels have therefore a better chance to contribute to 
the conductance in 3D than in Q1D. 

Figure presents the conductance distribution for 3D 
Anderson model at the critical point (W c = 16.5) and 
for 2D and Q1D systems with the same mean conduc- 
tance. As supposed, 3D distribution posses the longest 
tail, while two QlD-distributions are almost identical. 
From the plot of the -P(log g) (data not shown in the pa- 
per) we conclude that all four distributions are identical 
for small g. This is not surprise since the small -g be- 
havior is determined completely by the distribution of Z\ , 
which is the same in all systems. 

The quantitative difference between the conductance 
distribution for cubic and Q1D systems is measured by 
the integral 
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P(g)dg 



(5) 



In figure we plot 1(1) as a function of the mean conduc- 
tance (g) for various geometrical shapes of systems. For 
Q1D systems, we calculated both L x L z and L 2 x L z 
(L z » L). The data lay on the same curve. This is 
no surprise since the spectrum of z is linear in both sys- 
tems. A similar universality was found in 2D system. 
Here, we either increase disorder with constant system 
size or increase the size of square samples keeping dis- 
order constant. As supposed, the data for 3D and 2D 
systems scales independently 

The figure [| compares P(log<?) for 3D and Q1D sys- 
tems in localized regime. It confirms that P(\ogg) con- 
verges to Gaussian both for 3D and Q1D systems. The 
distribution for Q1D systems is much broader than that 
for 3D. While the shape of P(log<?) in weakly disor- 



dered Q1D systems can be obtained analyticallyOn-J, dis- 
tribution fo r . 3D systems is known only from numerical 
simulations£3l23 and deserves more detailed analysis. 



IV. STRONGLY LOCALIZED REGIME 

In the insulating regime all z increase linearly with the 
system length. The length-dependence of z\ determines 
the localization length £ 



, n 2L 
zi(L) = — + const 



(6) 



(For Q1D systems, L should be replaced by L z ). 

The knowledge of the spectrum of parameters z is 
important also in the insulating regime. As mentioned 
above, the spectrum-aLz remains linear in weakly disor- 
dered Q1D systemsQoO so that Z2 = 2z\ S> z\. Contrary 
to Q1D systems, we find that the difference (z% — Z\) in- 
creases only logarithmically with the system size for 2D 
systems (Figure 0): 



(z2 — z\) ~ c x logL 2D systems. 



(7) 



The same logarithmic behavior was found also in the 
2D systems with spin-orbit scattering (data are not pre- 
sented here). 

In 3D systems we find that the differences {z% — Z\) 
and (zs — zi) ■ . . do not depend on the system size in the 
strongly localized regime: 



(z2 — z\) w const 3D systems 



(8) 



In figure o we present the zi-dependence of the differ- 
ence (z2 — Z\) for the 2D and 3D systems. For z\ large 
enough, data scale to the same curve, confirming the one- 
parameter scaling. More important, data prove that the 
spectrum of the transfer matrix depends on the dimen- 
sion of the sample in the localized regime. 

Figure [7] presents the z\ - dependence of the differences 
(zi — Z\) for i = 2, 3,4 and 5. Although the convergence 
is worse for higher differences, the saturation to the Z\- 
independent constant is clearly observable as z\ increases. 

Due to (ft), the contribution of the second channel is 
nonzero whichever is the contribution of the first one. 
This contribution, however gives only small correction to 
the first channel contribution (logg) = log 4— (zi). Nev- 
ertheless, it is interesting to test how the contributions 
of the higher channels influence the statistical properties 
of the conductance. We found that apart the shift of the 
mean value the distribution P(— logg) and the distribu- 
tion P(zi) are almost identical. 

Although the absolute value of the conductance 
changes only slightly due to the contribution of higher 
channels, the present effect influences the variance of the 
log g. Indeed, the small difference z-i — Z\ causes that Z\ is 
confined from above: the probability to find z\ » (z\) 
is much less in the 3D than in the Q1D. The width of 



the distribution of z\ and of log g is therefore narrower 
in 3D than in in the Q1D. This agrees with numerical 
data presented in figure [|. Note also an asymmetry of 
the distribution: the probability to find log g — (log g) + 6 
is higher than the probability to find log g = (log g) — 8 
(<5>0). 

To support one-parameter scaling by more quantita- 
tive arguments we analyzed also statistical properties of 
the parameters z. Contrary to weakly disordered Q1D 
systems, we have to take into account also mutual corre- 
lations of Z\ and higher z in 3D. It is therefore not clear in 
advance, whether or not the one-parameter scaling the- 
ory is valid in this regime. However, Figure shows 
that the differences (z» — Z\) are unambiguous functions 
of (zi). Validity of the one-parameter scaling in the lo- 
calized regime is confirmed also by the data in figure o 
which shows that var z\ is an unambiguous function of 
the mean value (z±) in 2D even when disorder is strong. 
The same seems to hold also for 3D systems (figure g) 
In figure H we present also the (logg)-dependence of var 
log g. In the last case we find deviations from scaling 
only when the localization length £ becomes compara- 
ble with inter-atomic distance (£ rs 1). We analyzed 
also systems with non-zero Fermi energy. This was mo- 
tivated by Ref.E2l, where separate var z vs (z) - relations 
were found in systems with different Fermi energy. Our 
present data, however, scale to the same curve as that for 
the zero Fermi energy (Figure |8|). However, we cannot 
guarantee that more detailed statistical analysis will not 
show small deviations from the one-parameter scaling, 
similar to that observed in one-dimensional systems^. 



V. CONCLUSION 

We have shown numerically that the shape of the con- 
ductance distribution depends on the system dimension 
when the system leaves metallic regime. Qualitative ex- 
planation of this difference is based on the known spec- 
tra of the parameters z. In weakly disordered quasi-one 
dimensional systems the spectrum of z is linear indepen- 
dently on the dimension and the system length. How- 
ever, in the cubic d-dimensional samples the spectrum 
depends on the dimension of the system when the disor- 
der increases. 

We presented and compared the conductance distribu- 
tions for Q1D and 3D systems in the intermediate (crit- 
ical) regime. For both of them we found non-analytical 
behavior in g = 1 and exponential decrease of P(g) for 
g > 1. The distributions differ only quantitatively. 

We have shown for the first time that the differences 
(z2 — z\) increase logarithmically with the system size 
in two-dimensional strongly disordered systems. In three 
dimensional strongly disordered systems the differences 
{zi — z\) and (23 — z\) are constant and independent 
on the system size and on the-strength of the disorder. 
Is agrees with previous resultsLJ. The contributions to 



the conductance from the first, the second and higher 
channels have therefore the same system-size dependence 
~ exp —z\. Although the contribution of higher channels 
to the mean conductance is small, the shape of the spec- 
tra of z influences strongly the form of the distribution 
of the logg, which becomes narrower in 3D than in the 
weakly disordered quasi-one dimensional systems. The 
one-channel approximation is therefore not sufficient for 
the entire description of the insulating regime in three 
dimensions. 

Our numerical data confirm that Z2, zj, together with 
var Z\ are unambiguous functions of Z\ . We believe there- 
fore, that, in spite of the non-linearity of the spectrum 
of z, the random matrix theory can be generalized to 
the description of the non-metallic systemsE!rE2l. Such 
a generalization would approve applicability of the one- 
parameter scaling theory both to the critical and the 
strongly insulating regime. 

This work was supported by Slovak Grand Agency 
VEGA, Grant n. 2/7174/20. Numerical data were par- 
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Computer Center of the Slovak Academy of Sciences. 
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FIG. 1. Conductance distribution P(g) for cubic samples 
L x L x L with L — 8. Disorder increases from W — 10.6 
(metallic regime, (g) — 1.481) to W — 21.5 (insulator, 
(g) = 0.073). Critical disorder is W c — 16.5 and mean con- 
ductance at the critical point is {g) c ~ 0.28. For each value 
of the disorder, we plot P(g) (dashed line) and compare it 
with the conductance distribution for weakly disordered Q1D 
system which posses the same mean conductance (solid line) . 
Size of Q1D system is L 2 x L z with L = 8, disorder W — 4 
and length is tuned to get required (g) . Fixed boundary con- 
ditions were applied in x and y directions. 10 5 samples were 
collected for most statistical ensembles, (a) - metallic regime, 
P(g) is Gaussian for both cubic and Q1D geometry. Distribu- 
tions differ from each other only in their width: var ^3d f» 0.27 
and var <?qid ~ 0.12. This agrees with theoretical data (0.31 
and 0.133 for 3D and Q1D, respectivelyu) . Figures (b) - (e) 
present P(g) in the critical regime: (g) < 1. The main differ- 
ence between conductance distributions is observable in re- 
gion g > 1. (f) shows the conductance distributions in the 
insulating regime. Note that there are no Q1D sample with 
conductance g > 1 while P{g) for cubes still possesses tail up 
to g ps 1.5 
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FIG. 2. Critical conductance distribution for 3D Anderson 
model (dashed line) and conductance distribution for square 
and quasi-one dimensional systems Lx L z and L 2 x L z which 
have approximately the same mean conductance. As it is sup- 
posed, the quasi-one dimensional systems have the identical 
distributions. 
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FIG. 3. Integral 1(1) = J P(g)dg as a function of (g) for 
different geometrical shapes of samples. Dotted vertical line 
indicates mean critical conductance for 3D Anderson model. 
7(1) decreases more quickly for Q1D systems than for squares 
and cubes. For the quasi-one dimensional systems data for 
systems Lx L z and L 2 x L z scale to one curve. This is caused 
by of the universal (linear) spectrum of z. Data for squares 
and cubes scale differently due to the dimension dependence 
of spectrum of zs. 
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FIG. 4. (a) Distribution P(logg) for 3D and Q1D systems 
displayed in Figure hi (f ) . Note that although systems have the 
same (g), they differ in (logg). This indicates that mean con- 
ductance is not relevant parameter in the insulating regime. 
The width of the distribution is larger for Q1D system, (var 
log g ~ 6 for the Q1D system, and S3 4.33 for 3D system), (b) 
the same for 3D and Q1D systems with (log 5) ~ —15.3. Dif- 
ferent width of distribution is clearly visible: var log g « 25.66 
for the Q1D system, and only ~ 12.6 for 3D system. This is 
consistent with previous numerical results!! 
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FIG. 5. System-size dependence of the differences (22 — 21} 
and (23 — 21) for 2D orthogonal systems. Three different values 
of disorder were used. Solid lines are logarithmic fits. 2D 
symplectic systems exhibit the same logarithmic dependence. 
Inset: Plot of var z\ vs. {21}. 
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FIG. 6. (2i)-dependence of difference (22 — z\) for 2D (full 
symbols) and 3D (open symbols) disordered systems. Solid 
line is the theoretical prediction (22 — zi) = (21) for the Q1D 
weakly disordered systems. 
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FIG. 7. 21-dependence of the differences An{zi — 21) for 
i = 2,3,4 and 5. Disorder W = 28, 32, 38, 45, 54 and 60, 
system size L — 6, 8, ... 16. Data for W = 16.5 and W — 20 
(L = 8) are also presented to show how As As increase from 
their "critical" values (W = W c ). For large z\, they converges 
to constant, which do not depend neither on the system size 
nor on the disorder. 
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FIG. 8. var logg vs —logg for various W and L. Solid 
line is linear dependence var logg = —logg. Data lie on the 
same curve with exception of some points for L — 6 and 8, 
which represent systems with localization length smaller than 
inter-atomic distance. For completeness, we calculated also 
var log g for some systems with non-zero Fermi energy (= 4.4 
and 2.5) with L — 12, full diamonds. 



